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Abstract  ') 

! 

The  problem  of  ion  notion  and  laser  cooling  in  a  radio 
frequency  quadrupole  trap  is  considered.  It  is  shown  by  a  simple 
numerical  calculation  that  for  ion  motion  in  the  absence  of  laser 
cooling  the  effect  of  the  time  varying  potential  ,ir(R)cosn-t v 
can  be  approximated  by  a  time  independent  effective  potential 
whenever  the  frequency  of  oscillation  of  the  electric  field,  si‘ 
is  sufficiently  greater  than  a  resonant  frequency, 

S).\~  eVo  4a  ,  where  is  the  electronic  charge,  V0  is  the 
maximum  potential  applied  to  the  trap  electrodes,  and  tA  is  the 
mass  of  the  ion. 

Equations  for  the  Hamiltonian  matrix  elements  in  the 
presence  of  laser  cooling,  in  the  Lambe-Dicke  limit  are  derived 
and  studied  to  determine  the  cooling  transitions.  The  equations 
of  motion  for  the  density  matrix  in  the  Lambe-Dicke  limit  are 
derived  using  a  simple  representation  of  spontaneous  emission. 


QUANTUM  THEORY  OF  ION  MOTION  AND  LASER  COOLING 
IN  A  RADIO  FREQUENCY  QUADRUPOLE  TRAP 


I .  Introduction 

Purpose  and  Scope 

A  radio  frequency  quadrupole  trap  coupled  with  laser  cooling 
has  been  used  to  trap  a  single  ion  and  hold  it  in  a  volume  with 
dimensions  of  2/y.m  indefinitely  (16).  The  purpose  of  this  thesis 
is  to  develop  the  theory  of  particle  motion  in  the  trap  and  to 
determine  the  mechanisms  of  the  laser  cooling  process. 

Section  II  of  the  thesis  explores  the  problem  of  ion  motion 
in  the  absence  of  laser  cooling.  In  particular,  conditions  for 
stable  trapping  are  examined.  Two  approaches  to  this  problem  are 
taken.  In  the  first  approach  an  effective  potential  is  derived. 
The  effective  potential  solution  is  derived  starting  with  the 
actual  time  varying  potential  in  Schrod inger ' s  equation,  instead 
of  starting  with  the  classically  derived  effective  potential  as 
all  of  the  theoretical  developments  up  until  now  have  done  (4). 
The  second  approach  postulates  a  wave  packet  solution.  An 


equation  for  the  width  of  the  packet  is  derived  and  explored  both 
analytically  and  numerically. 

Section  III  of  the  thesis  analyzes  the  laser  cooling  process 
particularly  during  the  final  stages  when  the  ion  is  in  the 
Lambe-Dicke  limit,  that  is,  when  the  width  of  the  ionic  wave 
packet  is  small  compared  to  the  wavelength  of  the  of  the  cooling 
laser.  The  cooling  transitions  are  determined  by  evaluating  the 
Hamiltonian  matrix  elements.  The  equations  of  motion  for  the 
density  matrix  elements  are  derived. 


Background 


In  1980  Neuhauser,  et  al.  reported  trapping  a  single  barium 
ion  in  a  Paul  radio  frequency  quadrupole  trap  using  laser 
cooling,  and  holding  it  in  a  volume  with  dimensions  of  2yum 
(16:1137).  The  trapped,  cooled  ion  can  be  used  for  very  high 
resolution  spectroscopy.  Alone  in  the  trap,  its  spectral  lines 
will  not  be  broadened  by  collisions  and  if  it  is  restricted  to  a 
region  smaller  than  X/4  ,  the  Doppler  effect  will  be  suppressed 
(6:55)  • 

The  Paul  trap  consists  of  three  electrodes:  two  cap 
electrodes,  shaped  like  hyperbolas,  with  a  ring  electrode  between 
them.  An  oscillating  electric  field  is  applied  to  the 
electrodes.  The  trap  is  illustrated  in  Figure  1.  The  potential 


created  within  the  trap  space  is 


Wt,t)  '"it  (R)  COS  Jit 


(1) 

where  V  is  the  electric  potential,  \4  is  a  constant,  and  Tl 
is  the  frequency  of  oscillation  of  the  electric  field  (6:56-61). 
The  origin  of  coordinates  is  at  the  center  of  the  ring  electrode 
with  the  "2  -axis  running  between  the  two  cap  electrodes, 
perpendicular  to  the  plane  of  the  ring  electrode. 

A  Paul  trap,  which  has  an  oscillating  electric  field,  was 
used  instead  of  a  trap  with  a  static  electric  field,  such  as  a 
Penning  trap,  because  static  traps  must  use  an  electric  and  a 
magnetic  field.  The  magnetic  field  causes  a  Zeeman  effect  which 
is  undesirable  for  high  resolution  spectroscopy  (6:55). 

In  the  trap  used  by  Neuhauser,  et  al.,  the  spacing  between 
the  two  cap  electrodes  was  0.5  mm,  and  the  inner  diameter  of  the 
ring  electrode  was  0.75  mm.  The  frequency  of  oscillation  of  the 
electric  field,  il  ,  was  277  x  18  MHz  and  the  frequency  of 
oscillation  of  the  ion  in  the  well,  V  ,  was  277 x  2.4  MHz.  The 
barium  was  injected  into  the  trap  from  an  oven  and  ionized  by  an 
electron  beam.  The  ion  was  thermalized  by  allowing  helium  buffer 
gas  into  the  trap  at  a  background  pressure  of  <10-^  Torr 
(15:233)*  A  laser  tuned  to  the  Ba+  line  at  493*4  nn  was  shone  on 
the  trap,  and  when  a  barium  ion  was  present  in  the  trap  it  could 
be  observed  visually  by  fluorescence  (16:1139)* 

The  493*4  nn  line  for  Ba+  is  a  transition  between  the  ground 
6 2  S  i  / 2  state  and  an  excited  6^P1/2  state.  For  the  purposes  of 


cooling  the  ion,  the  laser  was  tuned  to  a  frequency  below  the 
resonant  frequency  of  the  transition  by  an  amount  .  One-third 
of  the  tine,  the  ion  spontaneously  decayed  to  a  52D1/2  netastable 
excited  state  instead  of  back  to  the  ground  state.  A  second 
laser,  tuned  to  the  5  ^  D  -] /2 -6^p  -j /2  transition  at  649.9  nm  kept  the 
ion  from  sitting  in  the  metastable  state.  An  energy  level 
diagram  is  shown  in  Figure  2.  Both  laser  beams  entered  the  trap 
between  the  ring  electrode  and  one  of  the  cap  electrodes,  at  an 
angle  of  45°  to  the  2.  -axis. 

This  experiment  generated  a  flurry  of  theortical  activity. 
Several  approaches  to  this  problem  can  be  taken:  systemising 
light  induced  forces  (3)>  (7);  investigating  the  density  operator 
of  a  particle  in  a  light  field  (2),  (11);  perturbation 
calculations  (9),  (  1 8 )  ;  and  quantum  statistical  treatments  (10). 
Reference  (18)  deals  with  the  cooling  of  both  free  and  bound 
atoms,  and  considers  particularly  the  case  of  a  harmonically 
bound  atom.  References  (9),  (10)  and  (11)  deal  specifically  with 
harmonically  trapped  ions. 

The  treatments  of  (9),  (10),  (11)  and  (18)  have  been  applied 
to  an  ion  in  a  radio  frequency  quadrupole  trap.  It  was  assumed 
that  this  was  justified  because  classically  a  particle  in  a 
rap  idly  oscillating  potential  field,  nr(F0  cos  at  ,  can  be 
considered  to  move  as  if  acted  on  by  a  time  independent  effective 
potential,  V'V»VrZ/'/lj  fA  _fl  (  12:90-93).  A  fully  quantum 
mechanical  justification  of  this  assumption  has  recently  been 


given  (4)  and  will  be  followed  in  Section  II  of  this  thesis  for 
ion  motion  in  the  absence  of  laser  cooling. 

Reference  (11)  is  the  most  detailed  analysis  of  the  final 
stages  of  the  laser  cooling  process  in  a  harmonic  trap,  and  the 
density  matrix  equations  derived  here  will  be  compared  to  the 
density  matrix  equations  used  there. 


II.  Ion  Motion  in  the  Absence  of  Laser  Coolin 


Effective  Potential  Solution 

Cook  et  al.  have  developed  the  theory  of  particle  motion  in 
a  rapidly  oscillating  field  in  Reference  (-4).  That  treatment  will 
be  followed  here  for  the  particular  case  of  an  ion  in  a  radio 
frequency  quadrupole  trap. 

In  the  absence  of  laser  cooling,  the  forces  acting  on  the 
ion  in  the  trap  are  due  entirely  to  the  rf  electric  field.  The 
electric  potential  for  this  field  is  given  by  equation  (1). 
Schrodinger ' s  equation  for  the  particle’s  wave  function  in  this 
case  is 


ifc  V2'V,-+eX'(R)cosilt  (2) 

To  understand  the  dominant  effect  of  the  potential,  consider  the 
solution  to  equation  (2)  if  only  the  potential  term  is  kept.  The 
solution  to  the  equation  is  then 

exp  {-i  IGr  e^fR)l4rR,o)6fR»(3, 

Therefore,  the  principal  effect  of  the  electric  field  is  an 
oscillating  phase  factor,  0  .  A  solution  of  tne  form 


‘VffyO-  <F(Ri't)9(R,t) 


(4) 


is  assumed.  The  advantage  of  switching  to  *p  is  that  ^  is  a 
slowly  varying  function  of  time  compared  to  the  rf  field  (4).  The 
equation  for  (P  is  then 


+  (xN-YM^sinMU  <P 


(5) 


Substituting  'a(t  -CoS  1/it 

for  s'ip2ilt  and  making  the 


rotating  wave  approximation  of  setting  the  rapidly  oscillating 
terms,  S»n Set  and  COS  Q.Q-L  ,  to  zero  yields 


^1?  =  +  <P 


(6) 


Equation  (6)  is  the  wave  equation  for  an  anisotropic  harmonic 

p  _  ■fi' eVo 

M  _/L 


oscillator  with  characteristic  frequency  V  ~  *  7-r -"V Y  in  the  Y 


and  Y  directions  and  twice  that  in  the  2  direction.  The 
total  wave  function  is  given  by 


n\  k^/t) 


Q(R,t) 


(7) 


where  Ll' ,  ^  is  an  appropriate  harmonic  oscillator  wave 


function . 


The  matrix  elements  of  the  Hamiltonian  in  equation  (2) 
between  the  wave  functions  of  equation  (7)  are 

H,o = 

=/dJR  v' & (■£ v*+  £'b(f()«sjit)'u„.9 

Applying  the  Laplacian  operator  yields 

Hn„'=/d3 R  -ih  hv 

4  evo  cosJlt  (X^V1-  22 


Substituting  y%(  i-cos2i£t)  for  s\n2Slt 

gives 

Hnn'  =  (dSR  Vnf^V 2  7Jn  (*V+  4?)  U»~] 

t[d3R"uJ  fix-  pun- 

-  V„>  (1 

MU1 

+  eV0  cosilt  (iH  Y  -XI1)  ^n'"] 

The  first  group  of  bracketed  terras  is  just  (^hn,Cn.  If  only 
motion  along  the  X  -axis  is  considered,  equation  (10)  becomes 


Hnn'  =  ^o'R+rn^^VoCosilt+^co^IIi) 

t  -lmL& 

+ ionta  Wfc’th  fev^o^^cosajx  * )  ni) 

+<Snn'420^&^oC°5i)i'Ka'CO53'at^ 

-  i  E  3;r,.ilir~] 

‘(VU'IVJM  7 

There  are  allowed  transitions  between  the  P  and  0+1  levels , 


with  transition  frequency  iV  .  The  terms  which  drive  these 
transitions  have  frequencies  H  and  2-fi  .  A  resonant 
f  requency ,  -file  is  defined  as  f£-/T  eVo/M  .  When  il=JU  o 
r^ii„/vn  ,  one  of  the  driving  frequencies  equals  the 
transition  frequency,  iv  ,  and  strong  transitions  are 
expected.  Actually,  it  is  shown  in  Reference  (4)  that  for  any 
positive  integer  ^  ,  ilo/'/F 

is  a  resonance  frequency  and 

strong  transitions  occur.  Thus,  stable  trapping  should  occur 
only  for  values  of  JX  larger  than  -fU.  For  the  parameters  of 
the  experiment  of  References  (16)  and  (15)>  -XL—  l»^-TLo  .  The 
trapping  stability  of  these  parameters  will  be  examined  for  the 
wave  packet  solution. 

The  results  so  far  derived  differ  from  the  effective 
potential  solution  assumed  in  Reference  (11)  in  two  respects. 


First,  the  phase  factor, 


,  does  not  appear  at  all  in  the 


solution  assumed  in  Reference  (11).  Second,  the  harmonic 
oscillator  wave  functions  of  Reference  (11)  are  considered  to  be 
stationary  states  in  the  absence  of  laser  cooling,  and, 


obviously,  the  wave  functions  of  equation  (7)  are  not.  The 
effect  of  the  phase  factor  on  the  laser  cooling  process  will  be 
discussed  in  Section  III. 


Wave  Packet  Solution 


For  one  dimensional  motion,  a  wave  packet  solution  to 
equation  (2)  is 


'tit,*}-  exp 


(  12) 


are  real  and 


where  IM  is  the  mass  of  the  ion,  q#)  and  <S(t) 
correspond  to  <x>  and  <P>  ,  and  O-'(i)  and  fa)  are 
complex.  The  importance  of  GJ  is  that  the  width  of  the  packet, 
w  ,  equals  fWM  (RaM  •  (3  represents  normalization  of  an 

phase.  Substituting  equation  (12)  into  equation  (2),  equating 
coefficients  of  powers  of  ")C  ,  and  assuming  that 


yields 


^  ~  -J.e.Vo  COSili  (13a) 

cosSlrt  (i3b) 

6  =  ki)Cu  +  ~2  cVo  Q2  Cos  Jit  (i3c) 


y.T“.T,-w7>w-V.'W! 


The  equations  for  and  (j  mean  that  the  particle  moves 
classically,  that  is  d<X>/cK  -  (  p>/m  and  d<p>/dt  =  <f>  , 

where  F  is  the  force. 

The  quantity  of  interest  is  G-’  .  The  equation  for  Ui  is  a 
generalized  Riccati  equation  (14:20-23).  This  equation  can  only 
be  solved  analytically  when  all  the  terms  are  polynomials.  If 
the  substitutions  U~-irn/fr)  and  T=  Jtt/l  are  made 
equation  (13b)  becomes 

dr^1  (4wlr)  cosiTl m  <  ’»> 


This  is  Mathieu's  equation.  Much  is  known  about  the  solutions  to 
Mathieu's  equation,  Mathieu  functions,  but  they  can  only  be 
represented  by  simple  functions  in  special  cases  (1:721-750).  The 
general  form  of  Mathieu's  equation  is 

'(r)U-  (a~^%  cosi'O'm  (  , 


For  the  values  of  and  A  used  in  the  experiment  of 

References  (15)  and  (  1 6 )  C\~  O  and  The  solution  of 

equation  (15)  for  these  values  of  and  is  in  what  is  known 

as  a  stable  region.  In  this  region  the  modulus  of  the  Mathieu 
function  is  always  finite  for  finite  times.  For  -Q-_T10, 
and  <^-2  |,4  I  •  The  solution  of  equation  (12)  for  these  values  is 


» wrjy 


*  ,'v.  vr:  = 


in  an  unstable  region.  In  this  region  the  modulus  of  the  Mathieu 
function  becomes  infinite,  on  a  finite  interval,  at  least  once 
(1:728).  If  ffl  ~  CX  +  |  b  then  the  width  of  the  wave  packet  is 

given  by 


/FfV+b*) 

W-  J  M  (a  V?'-  a'b 


The  modulus  of  the  Mathieu  function  is  /  a2  +  b1  ,  so  the 
numerator  of  equation  (16)  will  become  infinite  at  least  once  in 
a  finite  time  for  -n-Xl0.  This  does  not  necessarily  mean  that 
the  width  becomes  infinite,  which  corresponds  to  the  ion 
escaping,  because  the  denominator  might  also  become  infinite.  To 
understand  more  about  the  behavior  of  the  width  of  the  wave 
function  it  is  necessary  to  examine  equation  (13b)  numerically. 

Returning  to  equation  (13b),  if  - f~  }  »  T  =  V 't  , 

,  and  then 


ijiyz 

hco*ft+  Uf-y1) 


( 17a) 


(  17b) 


Equations  (17)  were  integrated  numerically  by  a  Runge-Kutta 
method  (8:212).  The  calculation  was  done  on  a  Columbia  Model  1600 
Personal  Computer  using  the  Basic  language.  The  text  of  the 


program  is  in  the  Appendix.  Three  sets  of  values  for  the 


constants  were  used;  the  values  from  the  experiment  of  References 
(  15)  and  (  16)  where  _/L  —  ,  values  for  and 

values  for  Jl--  JL . 

Figure  3  is  a  plot  of  the  width  of  the  wave  packet  versus 
tine  for  the  values  of  Vo  and  _Q  used  in  the  experiments  of 
References  (15)  and  (16).  The  initial  width  was  taken  to  be  the 


width  of  the  ground  state  harmonic  oscillator  wave  function, 
■w0=v/V'm7  .  The  width  oscillates.  The  oscillation  is 
actually  two  superimposed  oscillations  with  frequencies  Jl  and 
IV  .  The  width  stays  within  3  nm  of  the  original  value.  From  an 


initial  value  of  5.54  nm,  it  reaches  a  maximum  value  of  8.28  nm 


and  a  minimum  value  of  3*73  nm.  The  width  of  the  second  excited 


state  (first  allowed  transition)  of  the  harmonic  oscillator  wave 
function  is  given  by  N IB  W>.  For  this  case  that  would  be  12.4 
nm,  so  it  is  obvious  that  the  ion  is  not  being  excited  to  higher 


levels  and  is  stably  trapped. 

Figure  4  is  a  plot  the  width  of  the  wave  packet  versus  time 
for  the  value  of  V0  used  in  the  experiment  and  JL  —  lOjl^.  The 
initial  value  of  the  width  was  again  taken  to  be  the  width  of  the 
ground  state  harmonic  oscillator  wave  function,  in  this  case  12.6 
nm.  The  width  again  oscillates  with  two  frequencies,  JL  and 
IV  .  The  variation  from  the  initial  value  of  the  width  in  this 
case  is  much  less  than  the  variation  for  Jl=l,lJL  • 

Figure  5  is  a  plot  of  the  width  of  the  wave  packet  versus 
time  for  the  value  of  Vo  used  in  the  experiment  and  -11= Jl.. 


The  initial  value  of  the  width  was  again  taken  to  be  the  width  of 


the  ground  state  harmonic  oscillator  wave  function,  for  this  case 
3.98  nm.  The  minimunas  and  maxinums  for  this  case  occur  at 
regular  intervals.  The  spacing  from  minimum  to  minimum  and 
maximum  to  maximum  is  about  210  ns.  Each  maximum  is  about  145  ns 
from  the  minimum  preceeding  it  and  about  65  ns  from  the  minimum 
following  it.  There  is  a  maximum  of  1702  nm  at  611  ns,  which  is 
past  the  end  of  the  plot.  There  should  be  a  minimum  at  about  656 
ns.  The  value  at  this  minimum  could  not  be  calculated.  The 
program  calculates  the  values  of  y  and  from  equations  (14). 

The  width,  W  ,  is  \J  ^VjsA  y  »  so  that  a  minimum  in  W  is  a 
maximum  for  y  .  When  the  program  approaches  the  maximum  in 
at  656  ns,  the  value  calculated  for  y  overflows  its  register. 
This  corresponds  to  a  value  for  the  width  of  the  wave  packet  that 
is  less  than  10“20  cm.  This  is  an  unreasonably  small  value  for 
the  width,  so  the  overflow  is  probably  due  to  an  instability  in 
the  integration  technique.  However,  when  the  time  step  is  cut  in 
half,  the  value  for  y  overflows  at  the  same  point  in  time.  A 
further  examination  of  the  width  for  this  case  could  be  done, 


however,  it  seems  probable  that  the  particle  is  undergoing  strong 

transitions  for  the  chosen  values  of  V  and  n,  and  the  trend 

*  £ 

in  the  width  for  each  successive  maximum  to  be  larger  than  the 
one  preceeding  it  indicates  that  the  particle  will  escape. 


III.  The  Laser  Cooling  Process 


Hamiltonian  Matrix  Elements 


The  ion  in  the  trap  is  considered  to  be  a  two  level  system 
as  regards  electronic  transitions  for  this  section  of  the 
thesis.  The  basis  states  of  the  complete  system  are 


iin(R)  tir) 


rt«5£y* 


where  is  the  wave  function  for  the  electronic  state  with 

m-Gj  L  ,  is  the  wave  function  for  the  harmonic  oscillator 

state  with  fl  r  ^  t  *  /  <2. ;  ii>  ,  and  Q  is  the  phase  factor 

from  equation  (3).  The  energy  of  the  system  is  the  energy  of  the 
electronic  state,  plus  the  energy  of  the  harmonic 

oscillator  state,  fcV  (ft  4^)  .  The  energy  level  arrangement  is 

shown  in  Figure  6. 

The  Hamiltonian  for  an  ion  in  the  trap  in  the  presence  of 
laser  cooling  can  be  divided  into  two  parts,  the  Hamiltonian  for 

A 

the  internal  motion  of  the  outer  electron,  Hr  ,  and  the 
Hamiltonian  for  the  translational  motion  of  the  ion  in  the  trap, 

A 

Hr.  Each  of  these  can  in  turn  be  divided  into  two  parts,  the 

C  *  o  A  o 

part  that  is  present  in  the  absence  of  the  laser,  |~jr  and  Hr  , 

Ay  A  , 

and  the  part  due  to  the  laser,  H--  and  H-.  These  are 


[ij 


Energy  levels  and  transitions  of  the  ionic  system  in 
the  quadropole  trap.  The  electronic  level  is  indexed 
by  m,  and  the  harmonic  level  is  indexed  by  n. 


h;  =  -p(r')-£YR>i) 

H?  =  ^  4  aV, 

R*£-  -  ,^r  ((P^A(R,t)+MRt)-iP) 


(  19b) 
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where  p  is  the  momentum  operator  corresponding  to  the  outer 
electron's  position  with  respect  to  the  nucleus  of  the  ion,  V~  , 
fp  is  the  momentum  operator  corresponding  to  the  ion's  center 
of  mass  position  in  the  trap,  R  ,  rn  is  the  mass  of  the 
electron,  M  is  the  mass  of  the  ion,  El  is  the  electric  field 
due  to  the  laser,  and  fj  is  the  vector  potential  due  to  the 

A 

laser.  In  choosing  the  dipole  approximation  has  been  made. 

At  the  center  of  the  trap,  the  laser  radiation  can  be 
approximated  by  a  plane  wave,  so  that  £7  and  A  are 


F  =  E.coi(k'-R  -bjl) 
A  -  W  Fo  Sin  (k-  R  -  bit) 


The  matrix  elements  for  all  the  parts  of  the  Hamiltonian  of 
equations  (19)  with  the  wave  functions  of  equation  ( 1 8 )  can  then 

A 

be  found.  For  Hj  the  matrix  elements  are 

He  nmiVro  '  l  M x  I  Vnw)  '  £  XiP*™  (22) 
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For  Y\j  the  matrix  elements  are 

fy'nmn  W=-(li3r£pl?)  &H/d3R 


Three  simplifying  assumptions  can  be  made  about  ^  'f  13 
assumed  to  be  real  and  given  by  jj~  ^  0  where  X*  is  a  unit 

vector,  £  is  the  electronic  charge,  and  0.3  is  the  Bohr  radius; 
E0.p  -  Jj  5  and  a  two  by  two  matrix,  defined  below 

can  be  used  to  represent  the  integral  of  Wy  t  ,  over  . 


1  c>  1  \ 

l---; 

1  & 

f’mf o'  ~  (“*  ^ 

(.  o) 

(24) 
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Then  the  matrix  elements  of  Hi  can  be  written  as 

HjrnmnW  ~  -  EjwJ  d3R  U*  CO_s  (  fc.R-o^-U^ 

In  the  Lambe-Dicke  limit  is  much  less  than  unity  and 

Co>  s  ( lR.  R-  U'T )  -  c<os  u,+  4  V'  R  s\yn  cjt 

so  that  hjHWflW  becomes 


hi^nm n'rrt' ~  Eo inn '  COS  Lot 

^ i ooj-fc  h-  jdJR  i<Z  Run, 

For  one  dimensional  motion  of  the  center  of  mass  along  the 
^  -axis,  this  is 

^x^mn,rv)i  ~  Eoppim'  1 

-  M  tot 

X  nn'-| 

The  matrix  elements  of  Ji£  are 

nmriW  “  fjJR'U^fSva+'lr(R)c<«ni|<Uti,e 

Applying  the  Laplacian  operator  yields 

+  ^sinil£ 

-reV^co^iitCx^Yl^g1)  -yh.] 


Subst itut ing  '6  ( i  -  co  siOi)  for  ^,r>  2Jit 
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+  e.VaOasYL'fc  (X  +  Y-il  /"Uiv  J 


For  one  diraens 


■^coS2Ht) 


axis,  this  is 

H£  nmn'in'  =^m/  [in' fe4  ^2W\  (^CoS^4’^j^zc 

5'n-^1 
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It  was  shown  in  Section  II  of  this  thesis  that  in  stable  trapping 
situations  strong  transitions  between  the  H  andfflt^.  levels  do 
not  occur.  It  will  be  assumed  for  simplicity  that  no  transitions 
occur.  When  the  terms  involving  transitions  can  be  neglected, 
all  other  terms  involving  CO^-0-~fc  f  Cl  t  f  ancj 

coslilt  can  be  neglected  also,  so  that  H£nmhW  can  be 


•-4 


approximated  as 


V 

The  matrix  elements  for  are 

HfnmnW  -  S,w  Rl»n^ 

X  fv  R-wt)  Un<Q] 

+  -S>r>  (E'  R" W‘t')V('Un'Qi\ 

Applying  the  gradient  operator  and  noting  that  l 


&  E „=c 


becomes 


H^Wo-,'  '■Smw  (‘^)£0,Ij3R^  Sin(S.R-tJt) 

*£u„»  2 Off+Y-5-J2k)  ^  vo^„ri 

In  the  Lambe-Dicke  limit  .  R  is  much  less  than  unity 


that 


s  in  (  K-  R-^t)  =  cot  -  cot 


so  that  equation  (35)  becomes 


hEnmn'rfl'  = 

fj~  Sitiidt  jd^R  7^n("^r>'  ^ 

■fCo^Cot  ^r+V2>n/ 1] 


This  group  of  matrix  elements  cannot  be  as  simply  converted  to 
one  dimensional  motion  along  the  ^-axis,  but  this  is 
unnecessary.  The  interaction  between  the  ion  and  light  field 
described  by  this  part  of  the  Hamiltonian  is  nonresonant.  The 
cross-section  for  this  interaction  is  the  Thomas  scattering 
cross-sect  ion ,  °V=  eV^TT^M*^  .  For  the  intensity  of  light 

used  in  the  experiment  and  a  Ba+  ion  the  rate  for  this 
interaction  is  on  the  order  of  10“^  per  second,  which  much  less 
than  the  spontaneous  emission  rate  of  1.9  x  10^.  This  part  of  the 
Hamiltonian  can  be  neglected. 

It  can  be  seen  from  equations  (32)  and  (37)  that  the  phase 
factor,  0  ,  changes  the  matrix  elements  for  14"£  and  from 

what  they  would  be  if  Q  were  ignored.  However,  in  all  cases 

A 

\-j^.  can  be  neglected  entirely,  and  in  cases  of  stable  trapping 
the  terms  in  H £  nnnn'on1  caused  by  Q  can  be  neglected.  Thus, 
when  considering  laser  cooling  of  a  stably  trapped  ion,  the 
effect  of  the  phase  factor  can  be  neglected. 

The  matrix  elements  for  the  total  Hamiltonian, 

H-  B?  +  H|  +h;  are 

=  S»»m<Snn.rWm  +  *v(n*yJ] 

^osu't  (38) 
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■\>'  There  are  induced  transitions  between  the  THrO  ,  H  state  and 

the  (V)  —  ^  ,  D  or  D  t.  i.  state.  If  the  laser  is  tuned  below  the 
resonance  frequency  to  CJc~  \)  then  the  induced  transition 
between  the  ^  =  0,^  state  and  the  me  J.  ,  \1—  1  state  will 

dominate . 

Spontaneous  emission  has  not  been  included  in  the 
Hamiltonian.  The  dominant  spontaneous  transition  will  be  from  the 
¥>1—^  ,  n  level  to  the  fbrQ,  H  level.  Both  the  induced  and 
spontaneous  transitions  are  shown  on  the  energy  level  diagram  of 
Figure  6.  It  can  be  seen  from  the  diagram  that  the  net  effect  of 
the  combination  of  the  induced  transitions,  both  up  and  down,  and 
the  spontaneous  transition  down  will  be  a  decrease  in  the 
harmonic  oscillator  quantum  number,  which  corresponds  to 

CL 

'  cooling. 

It  should  be  noted  that  if  the  ion  is  cooled  to  the  ground 
state  of  the  harmonic  oscillator,  the  induced  transitions  will 
not  be  able  to  occur.  In  the  experiment  of  References  ( 15)  and 
(16),  the  ion  fluorescence  due  to  the  induced  absorption  and 
spontaneous  emission  can  be  observed  visually.  If  the  ion  is 
cooled  to  the  ground  state,  a  lessening  in  the  intensity  of  the 
light  should  be  observed.  The  light  will  not  disappear  entirely 
unless  the  natural  linewidth  of  the  transition  is  much  less  than 
the  harmonic  oscillator  energy  level  spacing.  The  natural 
linewidth  of  the  transition  used  in  the  experiment  is  277 x  19  MHz 
(5:11),  while  the  energy  spacing  between  the  harmonic  oscillator 
.  levels  is  27 7x  2 .  t}  MHz.  A  dip  in  the  fluorescence  could  be  used 


experimentally  to  determine  if  the  ion  has  been  cooled  to  the 
ground  state.  To  find  the  rate  of  cooling,  i.e.  the  rate  of 
change  of  ID,  and  to  determine  the  change  in  the  fluorescence 
when  the  particle  is  in  the  ground  state  it  is  necessary  to 
develop  the  density  matrix  equations. 


Development  of  Density  Matrix  Equations 


Z' 

A  density  operator,  ^  ,  can  be  defined  for  any  quantum 
system( 17 : 199) .  The  density  matrix  for  a  particular 
representation  is  ^  1  i  p  I  •  For  the  representation  used 

here  for  the  ion  system  this  is 


(39) 


n 


The  density  matrix  is  useful  to  work  with  at  this  point 
the  probability  Rim  that  the  system  is  in  a  state  ^ 
frimnrn’  If  the  density  matrix  is  known  at  any  time,  it 
found  at  later  times  by  the  relation 


because 

is 

can  be 


The  Hamiltonian  matrix  elements  should  be  for  the  complete 
Hamiltonian,  including  spontaneous  emission.  To  include 


spontaneous  emission  explicitly  in  the  Hamiltonian  it  is 
necessary  to  quantize  the  electromagnetic  field.  Spontaneous 
emission  can  be  included  in  the  density  matrix  equations  by 
following  the  treatment  of  Reference  (13).  The  density  matrix 
equations  including  spontaneous  emission  terms  are 

-for  m 4 m  '  or  n  ^  n ' 

^  (Horn  f  i  jio'ro'  -  n'wr) 


pn  m  n  ‘On  '  " 


'X  *+  r  ,  X  )  Pnmn 

■for  m r  m '  ahd  o 1 

*  _  l 

y  n  rn  n  m  "  )*£, 


ffY)  ‘ 


(41a 
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For  the  purposes  of  this  thesis,  spontaneous  emission  is  assumed 
to  occur  only  from  the  m  =  J,  ,  H  state  to  the  t  O  state. 

This  is  a  valid  assumption  if  the  natural  linewidth  of  the 
transition  is  much  less  than  the  spacing  of  the  harmonic 
oscillator  energy  levels.  For  this  case, 


^nmn'ro'  ^  Snn'  6  0 


(42) 


A  matrix,  n rv>  n ' m '  is  defined  as 


n  no  p  rn 


^hmin  p,  •  n.m#  "  I?,ynu  ^  i  -i  n'iv>'| 


f  4R1 


Substituting  the  matrix  elements  of  equation  (38)  into  equation 
(43)  and  noting  that 


?-  Shi-I  v/rT>  j\^n'nn  '  -  ^  (n+ ,)  _j  n '  jy,  / 

^in'4»  ]r>tn  -  sfr)'+  T 
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For  given  values  of  fn 
performed . 
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the  sums 
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can  be 
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Noting  that 
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the  density  matrix  equations  can  be  written  as 
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The  equations  for  the  density  matrix  elements  when  TO  TO  J 
contain  the  resonant  frequency  for  the  transition,  0^.,.  These 
elements  will  oscillate  rapidly.  A  substitution  can  be  made  to 
change  to  a  more  slowly  varying  function,  •  It  is 
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Jhe  equations  for  /  /  are 
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If  6iOCJ^  and  Ccx5tCt  are  written  in  terras  of  exponentials  it 


can  be  seen  that 
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because  c:  averages  to  zero  over  a  very  short  time  and  can 

be  disregarded.  For  cooling  purposes  the  laser  frequency , , 
is  tuned  below  the  resonant  frequency  of  the  transition,  by  an 
amount  U  .  So  that  the  equations  for  the  slowly  varying  parts  of 
the  density  matrix  are 
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These  equations  for  ^ncnn,m/  correspond  to  the  equations  for  the 
time  evolution  of  the  density  matrix  in  Reference  (11)  when  the 
assumptions  made  here  are  applied  there.  These  equations  of 
motion  for  the  density  matrix  elements  can  possibly  be  solved 
numerically  to  determine  the  cooling  rate  for  any  parameters,  and 
the  fluorescence  at  any  point  during  cooling. 


IV.  Conclusion 


In  summary,  the  problem  of  ion  motion  and  laser  cooling  in  a 
radio  frequency  quadrupole  trap  has  been  considered.  An 
effective  potential  solution  for  the  ion’s  wave  function  has  been 
derived  starting  with  the  full  time  varying  potential  in 
Schrodinger ' s  equation.  The  solutions  for  the  total  wave 
function  include  a  phase  factor,  0  ,  which  does  not  appear  in 
previous  solutions.  The  effect  of  this  phase  factor  is  to  cause 
transitions  between  the  IV  and  Hr  fv/+  2  harmonic 

oscillator  levels  of  the  particle  in  the  trap.  It  is  argued  that 
these  transitions  are  weak  and  can  be  neglected  for  values  of  the 

frequency  of  oscillation  of  the  electric  field,  -ft  sufficiently 
greater  than  the  resonant  frequency,  il0. 

A  wave  packet  solution  to  Schrodinger' s  equation  with  the 
time  varying  potential  is  presented.  An  equation  for  the  width 
of  the  wave  packet  is  derived.  It  is  shown  numerically  that  for 
A=i,9n.  the  width  of  the  packet  stays  close  to  the  initial 
value  of  the  width,  for  -fl  =  IQ-fl  o  the  width  of  the  packet 
oscillates  by  an  even  smaller  portion  of  the  initial  width,  and 
for  _fi=.Q_0  the  width  of  the  packet  increases  by  at  least  two 
orders  of  magnitude.  This  behavior  of  the  packet’s  width  was 
taken  as  confirmation  of  the  unimportance  of  transitions  for 
values  of  jft  greater  than  _fL  0  . 


The  Hamiltonian  matrix  elements  have  been  derived  for  a 


trapped  ion  with  two  electronic  levels  in  the  presence  of  laser 
cooling.  The  matrix  elements  are  derived  for  the  Lambe-Dicke 
limit.  When  the  laser  is  tuned  below  the  electronic  transition 
frequency,  L/J0  by  an  amount  d  ,  corresponding  to  the  harmonic 
oscillator  energy  spacing,  the  combination  of  induced  and 
spontaneous  transitions  causes  a  decrease  in  the  harmonic 
oscillator  quantum  number.  It  was  noted  that  when  the  ion 
reaches  the  ground  state  of  the  harmonic  oscillator  a  decrease  in 
the  observed  fluorescence  should  be  seen,  in  some  cases. 

The  equations  of  motion  for  the  density  matrix  have  been 
derived  using  the  Hamiltonian  matrix  elements  in  the  Lambe-Dicke 
limit.  These  equations  are  coupled  and  will  have  to  be  solved 
numerically.  Such  a  numerical  solution  is  planned,  and  the 
results,  if  available,  will  be  used  to  determine  the  rate  of 
cooling  in  the  Lambe-Dicke  limit  and  to  gauge  the  change  in 
fluorescent  intensity  when  the  ion  reaches  the  ground  state  of 
the  harmonic  oscillator  levels. 
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Appendix 


Computer  Program  for  Numerically 


Integrating  Equations  1 7 


10  '  Runge-Kutta  Calculation 
20  ' 

30  LET  DT= .0001 
40  LET  F=2 . 

50  LET  R=3.4344E-08 
60  LET  C=8 .2594E+07 
70  LET  Y=2.9202E+07 
80  LET  ZsO. 

90  LET  T=  0 . 

100  LPRINT  DT,  F,  R,  C 

110  LET  Wr.0021511*SQR( 1 . / Y) 

120  LPRINT  T,  Y,  Z,  W 
130  FOR  Ns  1  TO  120 
140  FOR  K= 1  TO  667 
150  LET  K 1 1  =  DT* ( 2 . * R* Y* Z ) 

160  LET  K12=DT*(C*C0S(F*T)+R*(Z~2-Y~2)) 

170  LET  K2 1  =  DT* (2.*R*(Y+K11/2.)*(Z+K12/2.)) 

180  LET  K22=  DT* ( C* COS (F*(T+DT/2. ) )  +R#  ( (Z+K12/2. ) ~2-( Y+K 1 1/2. )  ~2 ) ) 
190  LET  K3 1  =  DT* (2.*R*(Y+K21/2.)*(Z+K22/2.)) 

200  LET  K32=DT*(C»C0S(F«(T+DT/2.) )+R*( (Z+K22/2. )*2-( Y+K21/2. )*2) ) 
210  LET  K41=DT*(2.*R»(Y+K31)*(Z+K32) ) 

220  LET  K42=DT*(C*COS(F*(T+DT) )+R*( (Z+K32)*2-( Y+K31 )~2) ) 

230  LET  Y= Y+( 1 ./6. )*( K1 1+2 . »K2 1+2 . *K3 1+K4 1 ) 

240  LET  Z=Z+( 1 ./6.)*(K12+2.»K22+2.*K32+K42) 

250  LET  T=T+DT 
260  NEXT  K 

270  LET  W=. 002151 1*SQR( 1 ./Y) 

280  LPRINT  T,  Y,  Z,  W 
290  NEXT  N 
300  END 
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